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Abstract —Heterogeneous cloud radio access networks (H- 
CRANs) are potential solutions to improve both spectral and 
energy efficiencies by embedding cloud computing into heteroge¬ 
neous networks (HetNets). The interference among remote radio 
heads (RRHs) can be suppressed with centralized cooperative 
processing in the base band unit (BBU) pool, while the inter¬ 
tier interference between RRHs and macro base stations (MBSs) 
is still challenging in H-CRANs. In this paper, to mitigate this 
inter-tier interference, a contract-based interference coordination 
framework is proposed, where three scheduling schemes are 
involved, and the downlink transmission interval is divided into 
three phases accordingly. The core idea of the proposed frame¬ 
work is that the BBU pool covering all RRHs is selected as the 
principal that would offer a contract to the MBS, and the MBS as 
the agent decides whether to accept the contract or not according 
to an individual rational constraint. An optimal contract design 
that maximizes the rate-based utility is derived when perfect 
channel state information (CSI) is acquired at both principal and 
agent. Furthermore, contract optimization under the situation 
where only the partial CSI can be obtained from practical 
channel estimation is addressed as well. Monte Carlo simulations 
are provided to confirm the analysis, and simulation results 
show that the proposed framework can significantly increase the 
transmission data rates over baselines, thus demonstrating the 
effectiveness of the proposed contract-based solution. 

Index Terms —Heterogeneous cloud radio access networks, 
interference coordination, contract-based game theory. 

1. Introduction 

The demand for high-speed and high-quality data applica¬ 
tions is expected to increase explosively in the next generation 
cellular network, also called the fifth generation (5G), along 
with the rapid growth of capital expenditure and energy con¬ 
sumption m. Since the traditional third generation (3G) and 
fourth generation (4G) cellular networks, originally devised for 
enlarging the coverage areas and optimized for the homoge¬ 
nous traffic, are reaching their limits, heterogeneous networks 
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(HetNets) have attracted intense interest from both academia 
and industry 11. HetNets offer the advantages of serving 
dense customer populations in hot spots with low power nodes 
(LPNs) (e.g., pico base stations, femto base stations, small 
cell base stations, etc.), while providing ubiquitous coverage 
with macro base stations (MBSs), and reducing the energy 
consumption 0. Unfortunately, high densities of LPNs incur 
severe interference, which restricts performance gains and 
commercial applications of HetNets @1 . To suppress inter-LPN 
and inter-tier interference, coordinated multi-point (CoMP) 
transmission and reception has been presented as one of the 
most promising techniques in 4G systems O . Though perfor¬ 
mance gains of CoMP are significant under ideal assumptions, 
they degrade with increasing density of LPNs due to the non¬ 
ideal information exchange and cooperation among LPNs. On 
the one hand, it has been demonstrated that the performance of 
the downlink CoMP with zero-forcing beamforming (ZFBF) 
depends heavily on time delay, and thus the ZFBF design 
should be fairly conservative 0. Further, it has been shown 
that CoMP ZFBF has no throughput gain when the overhead 
channel delay is larger than 60% of the channel coherence 
time. On the other hand, it has been reported in ||2l that 
the average spectral efficiency (SE) gain of uplink CoMP in 
downtown Dresden field trials is only about 20% with non¬ 
ideal backhaul and distributed cooperation processing located 
on MBSs. Hence, novel system architectures and advanced 
signal processing techniques are needed to fully realize the 
potential gains of HetNets. 

Meanwhile, cloud computing technology has emerged as a 
promising solution for providing good performance in terms 
of both SE and energy efficiency (EE) across software de¬ 
fined wireless communication networks (H. By leveraging 
cloud computing technologies, the storage and computation 
originally provided in the physical layer can be migrated 
into the “cloud” to avoid redundant resource consumption 
and to achieve the overall optimization of resource allocation 
via centralized processing 0. As an application of cloud 
computing to radio access networks, the cloud radio access 
network (C-RAN) has been proposed to achieve large-scale 
cooperative processing gains, though the constrained fronthaul 
link between the remote radio head (RRH) and the baseband 
unit (BBU) pool presents a performance bottleneck ifTOll . Since 
the C-RAN is mainly utilized to provide high data rates in hot 
spots, real-time voice service and control signalling are not ef¬ 
ficiently supported. In order to avoid the significant signalling 
overhead through decoupling the user plane and control plane 
in RRHs, the traditional C-RAN must be enhanced and even 
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evolved. 

Motivated by the aforementioned challenges of both Het- 
Nets and C-RANs, an advanced solution known as the het¬ 
erogeneous cloud radio access network (H-CRAN) is intro¬ 
duced in mil- The motivation of H-CRANs is to avoid inter- 
LPN interference existing in HetNets through connecting to 
a “signal processing cloud” with high-speed optical fibers 
Ha. As such, the baseband datapath processing as well as 
the radio resource control for LPNs are moved to the BBU 
pool so as to take advantage of cloud computing capabilities, 
which can fully exploit the cooperation processing gains, lower 
operating expenses, improve SE performance, and decrease 
energy consumption of the wireless infrastructure. The simpli¬ 
fied LPNs are converted into RRHs, i.e., only the front radio 
frequency (RF) and simple symbol processing functionalities 
are configured in RRHs, while the other important baseband 
physical processing and procedures of the upper layers are 
jointly executed in the BBU pool. 

In H-CRANs, the control and user planes are decoupled, 
and the delivery of control and broadcast signalling is shifted 
from RRHs to MBSs, which alleviates the capacity and time 
delay constraints on the fronthaul. Therefore, RRHs are mainly 
used to provide high bit rates with high EE performances, 
while the MBS is deployed to guarantee seamless coverage and 
deliver the overall control signalling. With the help of MBSs, 
unnecessary handover and re-association can be avoided. The 
adaptive signaling/control mechanism between connection- 
oriented and connectionless modes is supported in H-CRANs, 
which can achieve significant overhead savings in the radio 
connection/release by moving away from a pure connection- 
oriented mechanism. 

As shown in Fig. [T] the BBU pool is interfaced to MBSs 
for coordinating the cross-tier interference between RRHs and 
MBSs. The data and control interfaces between the BBU 
pool and the MBS are SI and X2, respectively, which are 
inherited from definitions of the 3^^^ generation partnership 
project (3GPP) standards. With the help of centralized large- 
scale cooperative processing in the BBU pool, the intra¬ 
tier interference can be fully suppressed, while the inter-tier 
interference between RRHs and MBSs still present challenges 
to improving SE and EE performance. Significant attention 
has been paid to inter-tier interference collaboration ifTSll . ifT^ 
and radio resource cooperative management techniques ifTSll 
to alleviate inter-tier interference. As these works are focused 
primarily on HetNets, new models and methods are needed 
for H-CRANs. 

Intuitively, game theory is a useful tool for examining the 
inter-tier interference problem, and the published literature 
for HetNets has presented using this approach. For example, 
a Stackelberg game model has been used in uni to design 
a distributed radio resource algorithm, and its adaption to 
hierarchical problems such as hierarchical power control has 
been considered in |[T^ . More specifically, a Stackelberg game 
with a single leader and multiple followers is formulated 
to study joint objective optimization of the macrocell and 
femtocells subject to a maximum interference power constraint 
at the MBS in lfT9l . In II 20 II . a hierarchical game framework 
with multiple leaders and multiple followers is adopted to 



Fig. I. System model of the H-CRAN with one MBS and K RRHs. 


investigate uplink power allocation to mitigate inter-tier in¬ 
terference in two-tier femtocell networks. In order, for the 
Stackelberg game, to reach an equilibrium, players involved 
in the game have to interact periodically, which leads to high 
latency and signalling overhead. Furthermore, the incentive in 
the Stackelberg game model is ineffective if the information 
is asymmetric between the leader and follower. To overcome 
these issues, an advanced game model with proper incentive 
mechanisms and regarding the limited overhead is preferred. 

To overcome the aforementioned problems in the Stack¬ 
elberg model, contract-based game theory has recently been 
applied to aspects of spectrum sharing and relay selection. 
In ED, to solve the cooperative spectrum sharing problem 
under incomplete channel state information (CSI) in cognitive 
radio networks, a resource-exchange scheme is proposed and 
the corresponding optimal contract model is designed. The 
problem of spectrum trading with a single primary user selling 
its idle spectrum to multiple secondary users is investigated in 
EH, where a money-exchange based contract is offered by the 
primary user containing a set of different items each intended 
for a specific consumer type. In ||23l, under an asymmetric 
information scenario where the source is not well informed 
about the CSI of potential relays, a contract based game 
theory model is designed to help the source select relays for 
optimizing the throughput. Motivated by these existing works, 
since the cooperative processing capabilities and acquired CSI 
are significantly different between the BBU pool and the 
MBS in H-CRANs, contract-based game theory is a promising 
technique for modeling the inter-tier interference coordination 
mechanism under information asymmetry. 

Note that ideal inter-tier interference coordination depends 
critically on having accurate CSI for all transmitters and re¬ 
ceivers El; however, this assumption is non-trivial especially 
in view of the overwhelming signalling overhead and the 
inevitable existence of CSI distortion in practical systems ||25l. 
In previous studies, CSI errors are typically modeled as zero- 
mean complex Gaussian random variables without taking the 
practical training process and channel estimation methods into 
account ll^ . In order to present a realistic scenario, not only 
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the impact of imperfect CSI but also the corresponding training 
based channel estimation technique should be considered 
when examining contract-based interference coordination in 
H-CRANs. 

In this paper, we propose a contract-based interference 
coordination framework to mitigate the inter-tier interfer¬ 
ence between RRHs and MBSs in H-CRANs. Based on the 
proposed framework, different scheduling schemes can be 
utilized jointly to achieve performance gains. Specifically, 
three scheduling schemes are embedded into the proposed 
framework in this paper, and the downlink transmission is 
divided into three phases accordingly. At the first phase, 
RRHs help MBSs serve macro-cell user equipments (MUEs), 
which is termed the RRH-alone with UEs-all scheme. At the 
second phase, MBSs let RRHs use all radio resources to serve 
the associated UEs (denoted by RUEs), which is called the 
RRH-alone with RUEs-only scheme. Then at the third phase, 
RRHs and MBSs serve their own UEs separately, namely the 
RRH-MBS with UEs-separated scheme. Note that these three 
schemes are examples merely used to show the principles of 
the proposed contract-based inter-tier interference coordina¬ 
tion framework, and other inter-tier interference coordination 
schemes can be embedded into this framework, or replace 
these three schemes. Through the contract-based interference 
coordination framework, RRHs and MBSs establish a mutually 
beneficial relationship, which can be effectively modeled and 
optimized by using contract-based game theory. The main 
contributions of this paper are summarized as follows: 

• To coordinate the inter-LPN interference in HetNets, 
and to compensate for the substantial fronthaul overhead 
and processing in conventional C-RAN, the H-CRAN is 
considered in this paper, which incorporates advantages 
of both the HetNet and C-RAN. The H-CRAN decouples 
the control and user planes, and thus RRHs are mainly 
used to provide high data rates via cloud computing to 
suppress the intra-tier interference, while the MBS is 
deployed to guarantee seamless coverage and deliver all 
control and broadcast signals. 

• Since the cooperative processing capability and acquired 
CSI between the BBU pool and MBS are asymmetric, 
a contract-based interference coordination framework is 
proposed to coordinate the inter-tier interference between 
RRHs and MBSs, where the BBU pool and the MBS are 
modeled as the principal and agent, respectively. Based on 
the proposed coordination framework, three scheduling 
schemes, namely RRH-alone with UEs-all, RRH-alone 
with RUEs-only, and RRH-MBS with UEs-separated, are 
adaptively utilized. Accordingly, the downlink transmis¬ 
sion in each interval is divided into three phases. 

• With complete CSI, i.e., the BBU pool can acquire the 
perfect global CSI, sufficient and necessary conditions for 
a feasible contract are presented. An contract design that 
maximizes a rate-based utility is derived to achieve the 
optimized transmission duration for these three phases 
and to obtain the optimized received power allocation 
for RUEs and MUEs. Theoretical analysis indicates that 
the optimal rate-based utility is achieved when the time 


TABLE I 

Summary of Abbreviations 


AWGN 

additive white Gaussian noise 

BBU 

base band unit 

cdf 

cumulative density function 

CICF 

contract-based interference coordination framework 

CoMP 

coordinated multi-point 

C-RAN 

cloud radio access network 

CSI 

channel state information 

EE 

energy efficiency 

FRPC 

frequency reuse with power control 

H-CRAN 

heterogeneous cloud radio access network 

HetNet 

heterogeneous network 

IC 

incentive compatible 

IR 

individual rational 

KKT 

Karush-Kuhn-Tucker 

LPN 

low power node 

LSE 

least-squares estimate 

MBS 

macro base station 

MSE 

mean-square error 

MUE 

macro-cell user equipment 

OFDMA 

orthogonal frequency division multiple access 

pdf 

probability density function 

RB 

resource block 

RE 

radio frequency 

RRH 

remote radio head 

RUE 

remote radio head served user equipment 

SE 

spectral efficiency 

SNR 

signal-to-noise ratio 

TDD 

time division duplex 

TDIC 

time domain interference cancelation 

TTI 

time transmission interval 

UE 

user equipment 

ZFBF 

zero-forcing beamforming 


duration of the third phase is zero, i.e., the RRH-MBS with 
UEs-separated scheme is not triggered due to the large 
inter-tier interference even though the fairness power 
control algorithm is applied. 

• Considering the imperfect CSI acquired in piratical H- 
CRANs, training and channel estimation schemes are 
presented to obtain partial CSI. Based on the estimated 
partial CSI, contract design and optimization for the pre¬ 
sented framework are addressed. Specifically, sufficient 
and necessary conditions for a feasible contract under 
partial CSI are presented, which are individually rational 
and incentive compatible. In addition, the optimization 
problem to determine the contract that maximizes the 
rate-based utility is formulated. Moreover, an optimal 
contract design is derived when the number of RRHs is 
large. 

The rest of this paper is organized as follows: Section II 
describes the system model and the proposed interference 
coordination framework. In Section III, an optimal contract 
design with ideal and complete CSI is presented. Then in 
Section IV, contract design with practical channel estimation 
is considered. Numerical results are shown in Section V, 
followed by conclusions in Section VI. 

Notation: The transpose, inverse, and Hermitian of matrices 
are denoted by (•) , (•) , and (•) , respectively. || • || denotes 

the two-norm of vectors, diag (ei, 62 ,..., cat) is the x 
diagonal matrix. Ik represents the AT x AT dimensional unitary 
diagonal matrix. £{•} denotes the expectation of random 
variables and (•) denotes the first order derivative. Eor con¬ 
venience, the abbreviations are listed in Table HI 
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11. System Model and Intereerence Coordinated 
Framework 

The discussed H-CRAN consists of one BBU pool, K RRHs 
and one MBS, in which each RRH connects with the BBU 
pool through an ideal optical fiber, and the MBS connects 
with the BBU pool ideally with no constraints. It is assumed 
that all K RRHs cooperate, and they share the same radio 
resources with the MBS. The radio resources are allocated to 
different MUEs using orthogonal frequency division multiple 
access (OFDMA). It is assumed that each node has a single 
antenna, and thus for each resource block (RB) in the OFDM A 
based H-CRAN, only one MUE will be served by the MBS, 
and at most K RUEs will be associated with the K RRHs. 
Without loss of generality, M (M < K) RUEs and one MUE 
can share the same RB in OFDM A based H-CRAN systems. 
When M > K, K RUEs are selected to be auto-scheduled 
IIZTII . We note that a number of synchronization solutions 
for distributed antenna and virtual MIMO systems, e.g., the 
protocol in ll29l and post-facto synchronization (SOl, can be 
effectively applied to accomplish the synchronization among 
RRHs. Thus we assume that all RRHs and the MBS can be 
perfectly synchronized. 

Downlink transmission is considered, in which the BBU 
pool sends individual data fiows to RUEs via RRHs, and the 
MBS transmits the data fiows to MUEs. The radio channel 
matrix between RRHs and RUEs is denoted by G whose 
(m, k)-th entry gmk represents the channel coefficient between 
the k-ih RRH and the m-th RUE. The channel coefficient 
between the MBS and MUE is denoted by Let fkM denote 
the radio channel coefficient between the k-th RRH and MUE, 
and fBm represent the channel coefficient between the m-th 
RUE and MBS. It is assumed that all channels are quasi-static 
fiat fading so that they remain constant within one transmission 
block but may vary from one block to another. We further 
define Sm as the symbol transmitted from the BBU pool to 
the m-th RUE, and sm+i as the symbol from the MBS to 
MUE. It is assumed that distributed precoding is used at all 
K RRHs, and the K X M precoding matrix is denoted as A. 
The observations at all M RUEs can be expressed in vector 
form as 


yi? = GAs + fe-SM+l +11(7, (1) 

interference 

where s= [si, S 2 , • • •, sm]^, fs = [/bi,/b 2 , • • •,/bm]^, and 
11(7 is an M X 1 dimensional additive white Gaussian noise 
(AWGN) vector with each entry having the zero mean and 
variance cr^. Moreover, G = D 2 H models the independent 
fast fading and large-scale fading, where H is an M x AT 
matrix of fast fading coefficients whose entries are independent 
and identically distributed (i.i.d.) zero-mean complex Gaussian 
with unity variance, and D 2 is an M x M diagonal matrix 
with [D]^^ = Vm known at the BBU pool. The variance of the 
m-th element in is denoted by VBm- The signal received 
at the MUE is given by 

yB=gBSM+l-^ +^B, (2) 

interference 


where ^M = [/im, /2M, • • •, Ikm] with variances of elements 
denoted by um+i, and is zero-mean Gaussian noise with 
variance cr^. Eqs. o and ^ suggest that the received signals 
at the RUEs and MUE experience the inter-tier interference, 
thus degrading the transmission performance. 

To coordinate this inter-tier interference, three scheduling 
schemes are considered as examples, namely RRH-alone with 
UEs-all, RRH-alone with RUEs-only, and RRH-MBS with 
UEs-separated. Eor the RRH-alone with UEs-all scheme, both 
RUEs and MUEs are served by RRHs, and the MBS keeps 
silent. Eor the RRH-alone with RUEs-only scheme, to avoid the 
interference and interruption to RUEs, both MBS and MUE 
keep silent, and only RRHs and their corresponding serving 
RUEs operate. While for the RRH-MBS with UEs-separated 
scheme, both RRHs and the MBS operate simultaneously, and 
the interference between RUEs and the MUE is mitigated 
by using power control. To adaptively adopt these three 
schemes, a contract-based cooperative transmission framework 
is presented to coordinate these three scheduling schemes. 
Accordingly, as shown in Eig. [2l there are three phases for 
these three scheduling schemes in each time transmission 
interval (TTI) with a fixed time length Tq. 


- 

h 

To-ti-t2 

'-^ 



Phase I 

Phase II 

Phase III 



Fig. 2. Frame structure and scheduling algorithm 


• Phase I : During the time period of [0, ti] in each TTI, the 
RRH-alone with UEs-all scheme operates, where RRHs 
help the MBS serve the MUE, and the MBS benefits 
RRHs. Particularly, all K RRHs serve M RUEs and the 
single MUE in the typical RB, while the MBS remains 
idle. Since K > M -1- 1 is satisfied, the interference 
between RUEs and the MUE can be suppressed by using 
distributed precoding in the BBU pool. Note that RRHs 
consume additional power for serving the MUE, which 
can be regarded as the payoff of the RRHs. 

• Phase II : During the time period of [+, + +^ 2 ] in each 
TTI, the RRH-alone with RUEs-only scheme operates, 
where all K RRHs only serve the accessed M RUEs, 
while the MBS remains idle and the corresponding MUE 
is not served. Since RRHs help the MBS to serve the 
MUE in Phase I, the MBS awards RRHs by staying idle 
at this phase; thus there is no inter-tier interference from 
the MBS to RUEs. 

• Phase III : During the remaining time period of 
[U+^ 2 ,Lb] in each TTI, the RRH-MBS with UEs- 
separated scheme operates, where RRHs and the MBS 
serve RUEs and the MUE, respectively. To mitigate 
the inter-tier interference, the fairness power control is 
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utilized herein. 

The H-CRAN architecture can support efficiently the 
contract-based interference coordination framework, in which 
the BBU pool manages all connected RRHs and accessed 
RUEs. The inter-tier interference coordination happens be¬ 
tween the BBU pool and the MBS, which avoids the sub¬ 
stantial overhead and capacity requirements between LPNs 
and MBSs of HetNets. Note that the advanced inter-tier 
interference collaboration scheme ca and radio resource 
cooperative management techniques ca can be used to re¬ 
place the RRH-MBS with UEs-separated scheme at Phase 
III. To highlight characteristics of the proposed contract-based 
cooperative transmission framework and exploit performance 
gains resulting from contract-based game theory in H-CRANs, 
scheduling without coordinating the severe inter-tier interfer¬ 
ence is considered in the Phase III. 

Remark: The BBU pool and the MBS have conflicting 
objectives in the aforementioned three schemes. RRHs want 
to have a longer t 2 because all RUEs can be served without 
inter-tier interference during Phase II. Alternatively, the MBS 
wants a longer ti because the MUE can be served by RRHs 
without its own power consumption. If both RRHs and the 
MBS can benefit from this presented transmission framework, 
they would cooperate with each other. Otherwise, it would 
lead to ti = t 2 = 0, and ts = To which devolves to 
the traditional transmission scheme enduring severe inter-tier 
interference. The core idea of the proposed framework is that 
two participants under information asymmetry, i.e., the BBU 
pool and MBS, establish a mutually beneficial relationship, 
which can be effectively modeled and optimized by using 
contract-based game theory. 

III. Optimal Contract Design under Complete CSI 

It is preferred that BBU pool managing on all RRHs and 
RUEs act as the principal, due to the powerful capabilities 
of centralized cloud computing and large-scale cooperative 
processing. In the contract-based game model, the BBU pool 
as the principal designs a contract and offers it to the agent. 
The MBS as the agent decides to accept or reject the offered 
contract. Since the BBU pool and MBS are connected with 
an X2/S1 interface in H-CRANs, it is possible for them 
to exchange the corresponding individual CSI timely and 
completely. Thus, we can assume that the BBU pool can 
acquire the global and perfect CSI in this section. 

A. Rate-based Utility Definition 

A rated-based utility function for the contract is considered 
in this paper to represent the SE performance of H-CRANs. 
Though more advanced precoding techniques can be utilized 
to provide better interference coordination performance with 
higher implementing complexities than ZE does, the ZE pre¬ 
coding technique is adopted in this paper to suppress the inter- 
RRH interference in the BBU pool because this common linear 
precoding can achieve performance close to the sum capacity 
bound with a relatively low complexity when the number of 
RUEs is sufficiently large in interference-limited H-CRAN 
systems ca. 


• Phase I: Denote the received symbol power of Sm and 
<§M+i by Pci and Pmi, respectively. The sum transmis¬ 
sion data rate of all RUEs, denoted by Rci, and the data 
rate of the MUE, denoted by Rmi, can be expressed as 

Rci=M\og(^l + ^y (3) 

i?Mi=log(^l + ^), (4) 

respectively. Defining F= [G^, f^]^, the expected total 
transmit power of the RRHs is constrained by 

f {tr(FF^)“'A}<P,,ax, (5) 

where Pmax is the pre-defined allowable maximum trans¬ 
mit power, and A is an (M + 1) x (M + 1) diagonal 
matrix with the (M-hl)-th diagonal element given by Pmi 
and other diagonal elements given by Pci- The transmit 
power constraint can be re-written as 

tr(D-5£:{(HH")-i}D-U)<P,„ax, (6) 

where D = diag ... ,um+i) is an (M+1) x (M+1) 
diagonal matrix with H = [H^, hJ]^, where h/ is a 1 x P 
Gaussian vector with zero means and unity variances. As 
HH^ ^ W(M+1, Im+i) is an (M+1) x (M+1) central 
complex Wishart matrix with K{K > degrees of 

freedom, can be written as 

Y' _ ^^ _<p (7) 

• Phase II: Denote the received symbol power of Sm by 
Pc' 2 , and the sum rate of all RUEs can be expressed as 

Rc2 = Mlog(^l + ^y (8) 

where the long-term average transmit power is similarly 
constrained by 

"" (9) 

Since there is no inter-tier interference from the MBS to 
the desired RUEs due to the silent MUE, Pc 2 would be 
always set at the maximum value Pc 2 = to 

optimize the sum rate, where = X]m=i 

• Phase III: Denote the transmit power of the MBS and 
the RRH by Pb and Pcs, respectively. Under the fairness 
power control algorithm ||30l. Pcs and Pb are derived 
from 

{Pc 3 ,Ps} = arg max min {Pcs, Pms} , (10) 

Pc3,Pb 

s.t. 0<Pc3 ^ Pmax 5 0 <Pb ^ Pmax* ( 11 ) 
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The sum rate of all RUEs, denoted by and the data 
rate of the MUE, denoted by i?M 3 , can be expressed as 


Rm?> 


M 

Rcz = E] 

m=l 



C3 


l/Bmp-Pe + cr2 


( 12 ) 


log 


+ ^ 

Pc3fMG"(GG^)-'Gf^ +al)' 

(13) 


respectively. 

The rate-based utility of the BBU pool is defined as the sum 
rate in each TTI obtained by adopting the proposed framework 
compared with the traditional scheme in Phase III, which is 
given by 


Uc = tiRci + (^0 — h— ts) Rc 2 — {To — ts) Res- ( 14 ) 


TABLE II 


CSI 


Input; 


P^, where ^ [0, {K-M-l)vM+iPm 

satisfy > Pj]j, ?i = < 0, 


.(1) p(2) 






= P 


Ml 




Initialize: U^\ by substituting and P^\ into (14), respectively. 
k = 1, error e > 0. 


Repeat until convergence: 


Calculate: 


p(2) _p(l) 


Z = S - <il - <^ 2 , 


Pmi = P^m\ + (PmI - Pm\) (l 

If I < e 

Termination criterion is satisfied. 


W = \/ 

<;2+w+z A 

<?2—<?l+2ty ) ■ 


Eurther, the data rate-based utility of the MBS during each 
TTI can be expressed as 


Else 

Calculate Uc and = opmi ^Pmi=Pmi substituting Pmi into (14). 


Um = + hRMS- 


(15) If<^ = 0 

Termination criterion is satisfied. 


B. Contract Design under Perfect CSI 

To make the contract feasible under complete CSI, the 
principal should ensure that the agent obtains at least the 
reservation utility value that it would get by accepting this 
contract. Such a constraint is known as being individual 
rational (IR), which has the following definition. 

Definition {Individual Rational) A contract with the item 
(ti, ts, Pmi , Pci) is individual rational if the rate-based utility 
that the agent obtains is larger than the reservation utility u, 
i.e., 

tji^Ml + tsRMS >u = TqRmS- (16) 

Since the principal desires to obtain a maximal profit from 
the contract, the optimal contract (^i,^s,^ci) can be 
obtained from 

= max Urci + 

(To — ti— ts) Rc2 — (To — ts) Tesj, 

S.t. + I^sRmS > TqRmSi 

ts > 0, ti > 0, To — ti — ts > 0, 

I Pmi 

{K-M-l)ei {K-M-1 )vm+i~ 

(17) 

Lemma I: Under complete CSI, the contract 
{t\R%-,PMi,Pci) that maximizes the rate-based utility 
of the BBU pool satisfies the following conditions: 

tlRMl-{To-tl)RMS=^. (18) 

tl = 0. (19) 

p* p* 

_ 01 _^_ Ml _= p (20) 

Proof: See Appendix A. ■ 


Else 

If ^ < 0 

Pm\ = Pmi, = C and a = 
If ? > 0 

P^m\ = Pmi, = Pc and ?2 = 
k = k + l. 


Remark: Lemma 1 indicates that the MBS will obtain zero 
utility gain from the contract, which is not unexpected since 
the BBU pool will maximize its own utility with the minimal 
payoff. Moreover, the time period of Phase III would be zero 
when the contract is accepted by the MBS. Otherwise, the 
contract is rejected. 

Once the optimal contract is accepted, the objective function 
depends only on Pmi, and the computational complexity for 
obtaining P^^ is acceptable with numerical methods, which 
is provided in Algorithm 1. Once P^^ is obtained, t\ and P^i 
can be directly determined. 

Remark: Algorithm 1 is based on the spline interpolation 
method, whose computational complexity is mainly deter¬ 
mined by the number of iterations. In this case, the computa¬ 
tional complexity can be expressed as O (k), where k denotes 
the number of iterations. Note that k depends on the error 
convergence parameter e and the initial configuration of Pmi 
critically; thus it is important to choose proper initial Pmi to 
promote the convergence of the algorithm. 

Of course, complete and perfect CSI cannot usually be 
achieved in practice. In practical H-CRANs, CSI must be 
estimated and quantified, and is not error-free or delay-free. In 
the following section, we consider the impact of the imperfect 
CSI on the contract-based framework, where the principal can 
only acquire partial and non-ideal CSI via channel estimation. 
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IV. Contract Design under Practical Channel 
Estimation 


Similarly, the MSEs of and fuM^ denoted by and 
^/cM, can equivalently be written as 


In order to enable the BBU pool to acquire the downlink 
instantaneous CSI, an uplink training design is considered, 
which is based on the channel reciprocity in the time division 
duplex (TDD) mode. It is assumed that the fronthaul and 
backhaul are non-constrained and ideal, while the radio access 
links are non-ideal. 


A. Uplink Training Design and Channel Estimation 

Let denote the training sequence transmitted from the 
m-th RUE to the corresponding RRH, and V^m+i represents the 
training sequence from the MUE to the MBS. All have 
the N(> M) symbol length. ||'0^|p = A^Ps and ||'0 m+i|P = 
NPi), where Pg and Pi) are the training power of the RUEs and 
MUE, respectively. Moreover, ^ = [-01, -02, • • •, V^m] satisfies 
= NPs Im and = Om- During the training 

stage in the uplink, all RUEs and the MUE transmit their own 
training sequences simultaneously, and the observations at the 
k-th RRH can be expressed as 

M 

X/c Qmk'^m P f kM'4^M-\-l ? (21) 

m=l 


where is an N x 1 dimensional AWGN vector with 
covariance The least-squares estimate (LSE) of the 

channel coefficient between the k-th RRH and m-th RUE gmk^ 
denoted by Qmk^ can be expressed as 

and the LSE of the channel coefficient between k-\h RRH and 
MUE fkM is 


fkM 


V^M+l 

llV^M+lIP 


X/c. 


(23) 


Similarly, the MBS observes 


M 

Xb = '0M+1^H + + (24) 

m=l 


where wb is an N x 1 dimensional AWGN vector with 
covariance The LSEs of gs and /^m are 


9b = 


fBm — 


V^M+l 

IIV’M+lll 

lh^ 

'rm 


-xb, 


IIV’.i 


rXB, 


(25) 

(26) 


respectively. The mean-square errors (MSEs) of gmk and fBm^ 
denoted by 5mk and are equal and can be written as 


^mk — ^Bm ' 




IIV^M+lIP 

5i 

(27) 


SB=SkM = 


NP, 

S2 


(28) 


B. Contract Design under Imperfect CSI 

Eor notational simplicity, we define Rci, Pmi, Rc 2 ^ Rc^ 
and Rm under imperfect CSI to be compatible with the 
corresponding data rates Rci, Rmi, Rc 2 , Rc, and Rm under 
complete CSI. These transmission data rates under imperfect 
CSI are characterized in the following lemma. 

Lemma 2: The transmission data rates for all RUEs and the 
MUE at the three phases are given as follows. 

Phase I: The sum transmission data rate of all RUEs and 
the data rate of the MUE are given by 


Rci=M 

( 


log 


1 + 


Pcx 


M 

777^1 


H 


-1 


+ ^2Pm1 


FF 


H 


1-1 


-cr^ 


Jm+i,m+i / 

(29) 


Rmi - 


log 


1 +- 


Pmi 


M 

<5iPci^[fF 

777^1 


H 


1-1 


PS 2 PM 


FF 


H 


-1 




M+1,M+1 / 

(30) 


respectively, where the total transmit power is constrained by 


M 


y_i 

(K-M- 


Pci 


Ml 


m=l 


{K-M-l){v^ + Si) {K-M-1 ){vm+i + S 2 ) 


< P 

— max • 


(31) 


Phase II: The sum rate of all RUEs can be expressed as 

/ \ 


Rc 2 = M log 


1 + - 


C2 


V 


<5iPc2tr< 


GG^ 


1-1 


, (32) 
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where the total transmit power is constrained by 


M 

E 

777^1 


PC2 


{v^ + Si){K-M) 


<P 


Si -L max* 


(33) 


Note that RRHs prefer to achieve the maximal sum rate, 
and thus Pc 2 would always be set at the maximum value as 

With £2 = EEi 
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Phase III: The sum rate of all RUEs and the data rate of 
the MUE are given by 


Rc3 

M 

= 51 log 

m=l 


i+- 


PC3 


6iPc3tv< 


•| [gg^] ^ I 


> + \fBm\‘^PB + 0-l 


J 

(34) 


variable qi denoting the possibility of Obviously, 

qi e [0,1] and X] ® = 1- 

Hence, when the contract is accepted by the MBS with 
Igsf = 6’ the rate-based utility in the BBU pool can be 
written as 


Uc =t[Rci + (^0 - t[ - 4) Rc 2 

-{To-i) f p{z)Rc3dz. (39) 

J Z 


-RM3=l0g 


1 + - 


\9i 


?Pb 


PcsfG^ (gG") 


-2 


Gf"+(T2 


(35) 


Similarly, the rate-based utility in the MBS can be obtained 
by 

Um = t[RMi + 4^M3- (40) 


respectively. With the same power control algorithm as under 
perfect CSI, Pcs and Pb under the fairness power control can 
be achieved by 

{Pcs ,Pb} = arg max min | Res , Rms \ , 

Pc3,Pb t J 

S.t. 0 < Pcs < ^max, 0 < Pg ^ Pmax • (36) 

Proof: See Appendix B. ■ 

The BBU pool can acquire the estimated CSI of the RRH- 
RUE and RRH-MUE links, and it cannot estimate CSI of 
the MBS-MUE link, i.e., {gsf, which results in information 
asymmetry. In order to design the optimal contract, an in¬ 
centive compatible (IC) constraint is presented according to 
the revelation principle ED. It is assumed that Igsf has L 
potential quantified values and the set containing those values 
is denoted by E = {Ci, ^ 2 , • • •, Cu}- Without loss of generality, 
we further assume that < C 2 < • • • < • The IC constraint 

can be defined as follows. 

Definition {Incentive Compatible) A contract is incentive 
compatible if the agent with channel coefficients of the MBS- 
MUE link \gB^ = 6 prefers to choose the contract item 
(4,4? Tmi, Pci) designed specifically for its own type, i.e., 

^iPmI + 4^M3 — PmI + 4 ^M3 5 P} . 

(37) 

On the other hand, the IR constraint under partial CSI is 
given by the following definition. 

Definition {Individual Rational) A contract is individual ra¬ 
tional if the rate-based utility that the agent with = 6 

obtains by choosing the contract item (4:4? Tmi, Pci) is 
larger than the reservation utility ui , i.e., 

^iPmi + 4^M3 ^ — PqRms^'^^ ^ 2,..., L} , (38) 

where t{ and 4 are the time duration of Phase I and Phase 
III intended for the agent with {gB^ = respectively, and 
Rms calculated by (1351) under {gB^ = 6- 

Since the BBU pool can obtain the knowledge of the 
statistical parameters and the distributions of random variables 
by using appropriate learning and fitting methods, we assume 
that the BBU pool can acquire the probability density function 
(pdf) of 2 : = I/Bmp, the set S = and the 


Then the optimal contract with imperfect CSI can be ob¬ 
tained by solving the following optimization problem: 

max 

{{t{,tlPMi,Pci)yie{l,2,...,L}} 

Qi t[Rci + (To - 4 ~ 4) Pc2 - (To - 4) Pcs 

ZG{1,2,...,L} 

(41) 

s.t. 4^M1+4^M3 (42) 


t\RMl + t\pM3 — PmI + ^3 Pm3! 

Vl,l' 

4 >0,4 >o,To- 4-4 >o,w e {1,2,...,L}, 


Pmi > 0, Pci > 0, 
Pci 


(43) 

’(44) 

(45) 


Pmi 


{K-M-l)e 2 (K - M - 1 ) {vM+i + S 2 ) 


<P 


-1 -L max* 


(46) 


The above optimization problem can be further transformed 


to 


max Ucrr 

{(4,Pmi,Pci),V/ 6{1,2,...,B}} {(4,PMi.Pci).v!e{i.2,...,i}| 

(47) 


o<4 <4 < ■•• <4, 

(48) 

To-4-4>0,Wg{1,2,...,L}, 

(49) 

Pmi > 0, Pci > 0, 

(50) 

Pci Pmi 

< P_ 


{K-M-l)e 2 {K-M-1){vm+i+52) 


(51) 


where 

U*c 


{(4’^Ml>^Cl)>V(e{l,2, 


Pci — Pc 2 


^ qfii - R qit\ + ToR, 

ie{i,2,...,L} ie{i,2,...,L} 

(52) 


R = Rc2- Rc 3, i] = 

Rmi 


(53) 
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TABLE III 


CSI 


Input: where € 10,(K - M - 1)(vm+i + 

52)Pmaxl satisfy 

dU% I ^ , dU* I 

= SpfrlpM:=PU> < = ^lpMx=Pg> > 

Initialize: by substituting and into (|39), respec¬ 

tively, 

k = 1, error e > 0. 

Repeat until convergence: 

Calculate: 


3[[/*(2) _[/*(!)] _- ri) 

^ = -^(2)-fe-- ‘^2, W = \/- <^1<^2,Pm1 = Pmx + 


Ml Ml 

/p(2) _ p(l) \ / ^ _ S'2+'»^+^ ^ 

V^Ml ^2-^i+2'u;p 


If |p(2) _ p(l) I < 

Termination criterion is satisfied. 


Else 


Calculate Uc and c = I „ _ ^ by substituting Pmi into (|39j. 

'-^Ml—-P^Ml 


If ^ = 0 

Termination criterion is satisfied. 

Else 

If c < 0 

Pm\ = Afi, = Uc and ?i = ?, 

If ? > 0 

Pm\ = Pmi, = Uc and ?2 = ?, 

k = k + l. 


r — 


1 


UL 


Rmi 




(4 ^ ~ ^ 3 ) 


i=2 


phases under imperfect CSI can be rewritten as 

Pci 


(54) 


Note that the rate-based utility optimization in (l47l) with 
constraints dSli-dsB cannot always be solved in closed- 
form, which makes it difficult to design the optimal contract 
directly. Nevertheless, once the pdf of each random variable 
is determined, the optimal contract can be achieved according 
to the following Algorithm 2. 

Remark: The computational complexity of Algorithm 2 is 
similar to that of Algorithm 1. 

C Discussion with the Sufficiently Large K 

In this subsection, we consider the case in which the number 
of RRHs is large. It is assumed that each radio channel 
coefficient has a complex Gaussian distribution with mean 
zero. According to the Law of Large Numbers, we have the 
following approximations: 

LgqP 2FF^^^diag[D+,5ilM,(^^M+i+<52)] 

A A 

(55) 

Definition When the number of RRHs is sufficiently large, 
the transmission data rates of all RUEs and the MUE at three 


i?ci = Mlog 1 + 


SiP, 


ox. 


Ke2 


^2 Pm 


-RMl-log I 1+ 57pg— 

\ Ke2 


P('UM + l-h52) 

Pmi 


S 2 PM 1 


i?C2 = Ml0g(l + ^5^ 
\ Ks2 


M / 

Rc3 = J2^og 1 + ^ 

777^1 \ 


P(uM + l-h52) 

PC2 A 


+cri 


PC3 


Rm3 = 


= log 


1 + 


^ + \fBm\^PB+<jl 

B 


\9b\'^Pb 


VM + lPc3 




(56) 


(57) 


(58) 


(59) 


(60) 


V 

Eqs. ([56b-([6Qb in Definition 4 can be derived directly by 
substituting (1551) into the transmission data rates in Lemma 
2. Based on Definition 4, the optimal contract design can be 
derived via the following lemma when the number of RRHs 
is sufficiently large. 

Lemma 3: Eor a sufficiently large K under imperfect CSI, 
once the contract is accepted by the agent, the optimal fg, P^ 


Cl’ 


and Pmi follow 


4 = 0,/ = 1,...,L-1, 


( 61 ) 


^3 =< 


To 

[0,To] 

0 


Rmi{Rc 2 — Res) < RM 3 iRc 2 — Rci), 
Rmi{Rc 2 — Res) = Rm 3 {Rc 2 — Rci), 
Rmi{.Rg 2 — Res) > Rj(^^{Rc 2 — Rci), 

(62) 


p* 

^C1 




p* 

^M1 


= Pr. 


(63) 


{K-M-1)s2 {K-M-1){vmxiP^2) 

Proof: See Appendix C. ■ 

V. Numerical Results 

In this section, Monte Carlo simulations are provided to 
evaluate the performances of the proposed contract-based 
interference coordination framework in H-CRANs numeri¬ 
cally. All the fast fading channel coefficients are generated 
as independent complex Gaussian random variables with zero 
means and unit variances. The large-scale fading coefficients 
Vmim = l,...,{M+ 1)) are modeled as ■Um =2m/(?’m/'ro)", 
where Zm is log-normal random variable with standard devia¬ 
tion SdB, Vm is the distance between the m-th UE and RRHs, 
ro is set to be 100 meters, and 1 ; = 3.8 is the path loss exponent. 
The other common parameters are set at Pmax = 1, To = 1, 
M = 4 and = 10; thus the signal-to-noise ratio (SNR) is 
given by SNR= Totally 10^ Monte-Carlo runs are 

used for the average calcufations. Two baselines are considered 
in this paper: 

• Frequency reuse with power control (ERPC): RRHs and 
the MBS reuse the same RB, and the transmit powers 
are determined according to the fairness power control 
adopted for the proposed CICE in Phase III: 

{-Pc 3 ,-PB} = arg max min{Rc3,RM3], 

Pc3,Pb 

s.t. 0<Pc3 ^ Tmax: 0 <Pb ^ Pmax* 
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• Time domain interference cancelation (TDIC): RRHs and 
the MBS transmit separately in the different time phases, 
where the transmission duration allocated to RRHs and 
the MBS are Ir and Im, respectively, and they satisfy 
the relationship with tR = tM = ^- 

A. Performance Evaluations under Perfect CSI 

In Fig. [3l the sum rates of all RUEs versus SNRs under 
perfect CSI are compared among the CICF, FRPC and TDIC. 
According to these simulation curves, the proposed CICF 
achieves a larger sum rate than the two baselines do, which 
indicates that the BBU pool as the principal can benefit 
from the proposed contact model, thus demonstrating the 
effectiveness of the proposed CICF for all RUEs. Due to gains 
from the optimized time duration of the three phases and the 
allocated transmission power of RRHs and the MBS, the sum 
rate of all RUEs increases almost linearly with SNRs, whereas 
the baseline ERPC cannot provide the same performance gain 
as the CICE because the intr-tier interference is still severe 
even though the fairness power control is utilized. Eor the 
baseline TDIC, the sum rate is the lowest because its spectral 
efficiency is low even though the inter-tier interference is 
avoided. 



Fig. 3. Sum rates versus SNRs for all RUEs with different schemes. 

To evaluate performance gains of the agent, the data rates of 
the MUE versus SNRs are compared among the CICE, ERPC 
and TDIC in Pig. [H It can be observed that the proposed 
CICE achieves the equivalent data rate for the MUE as the 
ERPC, which is not unexpected since the equality of the 
IR constraint would always hold when the BBU pool can 
acquire the complete and perfect CSI. Comparing with the 
baseline TDIC, the data-rate based utility gain for the MUE is 
significant for the proposed CICE because only partial spectral 
resources are utilized in the TDIC. 

Summarizing the sum rates of all RUEs and MUE in Pig. 
[3] and Pig. [H respectively, the rate-based utility gain from 
the proposed CICE is significant, which demonstrating the 
effectiveness of the proposed contract-based solution under 
perfect CSI. 



Fig. 4. Data rates versus SNRs for MUEs with different schemes. 

B. Performance Evaluations under Partial CSI 

In Pig. m the sum rates of all RUEs versus SNRs under 
imperfect CSI for the proposed CICE and baselines (i.e., ERPC 
and TDIC) are evaluated and compared. Similarly to the results 
under perfect CSI, the proposed CICE can achieve a significant 
sum rate gain compared with the baselines due to the fact that 
both the time durations of the three phases and the allocated 
power are optimized. Moreover, comparing the performance 
gain gap from the CICE to ERPC in Pig. [3] and Pig. [5l the 
benefit of the principal under imperfect CSI is lower than 
that under perfect CSI. This happens because the derived 
optimal contract under imperfect CSI cannot always ensure the 
equality of all IR constraints. Purthermore, the performances 
of ERPC and TDIC under perfect CSI is better than that under 
imperfect CSI due to the fact that the channel estimation error 
degrades the precoding performance at RRHs. 



Fig. 5. Sum rates versus SNRs for RUEs with different schemes. 

In Pig. [6l data rates of the MUE versus SNRs under im¬ 
perfect CSI are compared among the CICE, ERPC and TDIC. 
Comparing with ERPC, the agent can achieve a significant 
performance gain in the proposed CICE because the optimal 
contract under imperfect CSI cannot ensure the equality of all 
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IR constraints. Meanwhile, the performance of FRPC in Fig.[6| 
is worse than that in Fig. [4| because the fairness power control 
depends critically on the instantaneous CSI. While TDIC in 
Fig.[6]is the same as that in Fig. [4] since the CSI has no impact 
on the downlink transmission of the MBS with this baseline. 



VI. Conclusions 

In this paper, a contract-based interference coordination 
framework has been proposed to mitigate inter-tier interference 
in H-CRANs, whose core idea is that the BBU pool is selected 
as the principal that would present a contract to the MBS, 
and the MBS as the agent decides whether to accept the 
contract or not according to its potential for improving the 
SE performance. An optimal contract design that maximizes 
the overall rate-based utility has been derived when perfect 
CSI is available at both principal and agent. Furthermore, the 
contract optimization under the situation where the partial 
CSI has been obtained from channel estimation has been 
considered as well. Monte Carlo simulations have been pro¬ 
vided to corroborate the analysis, and simulation results have 
shown that the proposed framework can significantly increase 
the sum data rates, thus demonstrating the effectiveness of 
the proposed contract-based framework. Interesting topics for 
further research include, more advanced inter-tier interference 
coordination schemes to enhance the RRH-MBS with UEs- 
separated phase in the presented framework, and the inclusion 
of fronthaul constraints within this framework. 


Fig. 6. Data rates versus SNRs for MUEs with different schemes. 

To characterize the optimal contract in the proposed CICF, 
the cumulative density function (cdf) of both t\ and are 
plotted in Fig. [T] It can be seen that P{t\ >0.5To) >0.7 and 
P(t2 <0.5To) >0.8, which indicates that the time duration of 
Phase I should be longer than that of Phase II to optimize the 
sum date rates. To optimize the transmission data rates of the 
whole H-CRAN, more time resources should be assigned to 
allow MUEs to be served by RRHs, and thus an example 
of the transmission frame structure under these scheduling 
algorithms and configurations in this paper is suggested such 
that the time duration of Phase I is about 70 percent of the 
TTI length, and the time duration of Phase II is about 30 
percent of the TTI length. Note that the transmission duration 
of Phase III would be zero when the contract is signed, which 
is indicated in both Lemma 1 and Lemma 2. 



Fig. 7. The cdf distributions of and 


Appendix A 
Proof of Lemma 1 

Lor the fixed ti and ^ 2 , the optimal Pci and Pmi can be 
derived by using the Karush-Kuhn-Tucker (KKT) conditions. 
The corresponding Lagrangian function can be written as 


L{Pci , Pmi , A) = -tiRci 
Pci 


-A 


P — - 

-L max 


Ml 


{K-M-l)si {K-M-1)vm+i\ ’ 
whose solution is derived as 

1 


(64) 


WL{Pci,Pmi,X)= -h^^+X 


Pci 


dPci ' "{K-M-l)ei 
Pmi 


= 0, (65) 


{K-M-l)ei {K-M-1)vm+i, 


A > 0. 


= 0 , ( 66 ) 

(67) 


It is obvious that A must be positive to achieve the optimal 
solution for Pci, and thus the optimal received powers P^-^ 
and allocated to RUEs and MUE should satisfy 

P* p* 

_ ^C1 _^_ -^Ml __ p f 

^K-M-l)ei {K-M-1 )vm+i ^ ^ 

Lor given Pci and Pmi, the optimal ti and t 2 can be 


obtained by solving following equations: 

Rc 2 -RciPv-vi- wRmi = 0, (69) 

i?C'2—+ =0, (70) 

Viti=0, V 2 t 2 = 0, v{T — ti—t 2 ) = 0, (71) 

Tu{tiRMi PtsRMs-u) = 0, (72) 

>0, V2>0, v>0, >0. (73) 


where v, vu, vi, and V 2 are the Lagrangian multipliers. It 
is obvious that m must be larger than zero. Therefore, the 
IR constraint is derived. Moreover, in order to maximize the 
transmission bit rate of the MUE, the optimal tg should be 
equal to zero. 
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Appendix B 
Proof of Lemma 2 

Denoting the m-th column of and A by and a^, 
respectively, the transmission data rate can be derived as ll^ 

logfl +__"i . 

\ cr2+var{g,„a„^} + |/Bm|2PB + E^rT^^[|g"^»i^/ 

(74) 

For the ZF precoder, we have gm^i = Smi- Denote the radio 
channel matrix by ^ = G+G, where G is the estimation error, 
and we have 




(75) 


According to (EHi, the expectation and variance can be 
written as 


M 


£[gmeii] = Pci, £[\gmeii\^] = S,Pci ^ [GG^],. , 

i^m 


var{gmai}=£ [a^g^gmam] =(5iPci [GG^ 


(76) 

(77) 


Thus, the sum data rate of all RUEs in the downlink 
transmission during Phase III is given by 

/ \ 

PC3 


M 


-Rc3 = L]l0g 


m=l 


i+- 




()iPc3tr< 


■| [gg^^] ' I 




J 


(78) 


and the data rate of the MUE can be written as 
\9b\^Pb 


Rm3 = log I 1 + 


cr2 + |fAr 




= log 


1 +- 


\9b\^Pb 


^ Pc3fG"(GG^^) ^Gf^ + 


(79) 


Similarly, the sum rates during the Phase I and II can be 
calculated in the same way. 


Appendix C 
Proof of Lemma 3 

By using the KKT conditions, the optimal can be 
obtained by solving the following equations: 

- {Rc2-Rci)QL 

TlMl 

+ (-^C2—-^C3)^L—^l+^2(l— 5^^ ) = O 5 (80) 

TlMl 

cjifg =0, uj2{To-i^-t^) = 0, (81) 

cji>0, CJ2>0. (82) 

Based on (ISOl) . (ISTI) . and (1821) . we have the following results: 

• When ^mi(^C 2 - Res) > ^m 3 (^C 2 must 
be positive, i.e., = 0. As tg (/ = 1,..., L — 1) is often 

less than tg, 4 is equal to zero. 


• When ^mi(^C 2 — Rc 3 ) = Rms{Rc 2 —Rci), ^3 can be 
an arbitrary value in the interval [0,To]. The derivative 
of —U^ with respect to can be written as 

{Rc2-Rci) [QL^QI^i] {Rm3 ~ ^M3^) ^ (^^) 

To maximize we have =0, and all (I < L — 1) 
should be equal to zero. 

• When ^mi(^C 2 - Rc 3 ) < ^m 3 (^C 2^2 must 
be positive. Meanwhile, the following equation holds: 

To=i^+tl (84) 

Obviously, this contract is unacceptable for the agent. 
Therefore, should be equal to Tq. In order to satisfy 
the IR constraint, tg (/ = 1,..., L — 1) should be set as 
zero. 

As for the optimal powers assigned to RRHs and MBS in 
Phase I, when 4 (/ = 1,..., L — 1) is equal to zero, can 
be expressed as 

^^^iR£^z3£3hlL+To{Rc3-Rc2). (85) 


The Lagrangian function can be written as 

{Rc1—Rc2)uL 


Li{Pci, Pmi, = - 


Rmi 


fl{Pci,PMl) 


. \p Pci Pmi 

1 ^ max {K-M-1){vm+1 + S2) ' 

( 86 ) 

The derivative of the Lagrangian function with respect to 

Pci is 


VI/i(Pci, Ai) = 

dfijPci^PMi) . f 

dPci \K-M- 1)62 ’ 

^ p_ Pg _ Pmi _ 

{K-M-1){vm+1 + 52) 

Ai > 0. 


(87) 

= 0 . 

( 88 ) 

(89) 


Since the derivative of fi{Pcc Pmi) with respect to Pci 
is positive. A* must be positive, and the optimal powers P^i 
and P^^ should satisfy 


Cl 


Ml 


{K-M-1)s2 (P-M-1)(um+i+^2) 


= Prr 


(90) 


If tg is positive, can be re-written as 


U^ = {Rci-Rc 2 )[l-qL] 
P{Rci—Rc2)qL 


Ul 
Rmi 
^L—Rm3 


Rmi 


pTo{Rc3-Rc2)[f - Ql]‘ 

(91) 






























IEEE JOURNAL ON SELECTED AREAS IN COMMUNICATIONS, VOL. 33, NO. X, SECOND QUARTER, 2015 


13 


Similarly, the Lagrangian function is re-formulated as 

L2{Pc1,Pm1,^2) = 

- (^C1-^C2)[1 - Ql] P{Rci-Rc2)qL 

Rmi Rmi 


— Ao 


/2 {Pci ,Pmi) 

Pci Pmi 


(92) 


{K-M-l)e 2 (i^-M-l)(uM+i + ^2) J 


(93) 


As the derivative of f 2 {Pci^ Pmi) with respect to Pci is 
also positive, the optimal powers allocated to P^i and 
are the same as in (l90l) . 
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